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Introduction 

In a number of studies it has been demonstrated that systems of the form 

{doM + A)u = F, (0.1) 

where A4 is a continuous, hnear mapping and the densely defined, closed linear operator 
A is such that A and A* are maximal w-accretive for some suitable a; G M, cover numerous 
models from mathematical physics. Indeed, A skew-selfadjoint0 is a standard situation, 
which for simplicity we shall assume throughout. The well-posedness of (10. ip hinges on a 
positive definiteness assumption imposed on OqM in a suitable space-time Hilbert space 
setting. Under this assumption the solution theory is comparatively elementary since 
{doAi + A) together with its adjoint are positive definite yielding that {doAi + A) has 
dense range and a continuous inverse. 

In applications of this setting the operator A has a rather simple structure whereas the 
complexity of the physical system is encoded in the "material law" operator Ai. A simple 
but important case is given by 

M = Mo + do^Mi, 

where Mq, Mi are time-independent continuous linear operators. Here we have antici- 
pated that in the Hilbert space setting to be constructed Oq^ (time integration) has a 



^Two densely defined linear operators A,B are skew-adjoint (to each other) ii A = —B*. li: A = B, we 
call A skew-selfadjoint (rather than self-skew- adjoint). The proper implications 

A selfadjoint A symmetric A Hermitean 

are paralleled by 

A skew-selfadjoint A skew-symmetric A skew-Hermitean. 

Frequently, "skew-adjoint" is used to mean "skew-selfadjoint". We shall, however, not follow this custom 
for the obvious reason. 



4 



proper meaning. The positive definiteness assumption requires Mq to be non-negative and 
self adjoint and 

gMo + $He Ml > cq (0.2) 

for some Cq G ]0, oo[ and all sufficiently large g G ]0, oo[. Since we do not assume that Mq 
is always strictly positive, (IU.2p may imply constraints on Mi. If Mq is positive definite it 
may seem natural, following the proven idea of first finding a fundamental solution (given 
by an associated semi-group), and then to obtain general solutions as convolutions with the 
data (Duhamel's principle, variation of constants formula) and so proving well-posedness. 
This is the classical method of choice in a Banach space setting, see e.g. [H HI |5l [H [T7] as 
general references. In comparison our approach is (currently) limited to a Hilbert space 
setting, however, apart from being conceptually more elementary, it allows to incorporate 
delay and convolution integral terms by a simple perturbation argument and, if Mq has a 
non-trivial kernel, the system becomes a differential-algebraic systems, which to the above 
approach makes no difference, but cannot be conveniently analyzed within the framework 
of semi-group theory. 

The purpose of this paper is to extend well-posedness results previously obtained for time- 
shift invariant material operators Ai to cases, where Ai is not time-shift invariant. This is 
the so-called time-dependent or non-autonomous case. The above-mentioned limitations 
of the semi-group approach carry over to the application of classical strategies based on 
evolution families intoduced by Kato ([7|), for a survey see e.g. [TB| [S], which are the 
corresponding abstract Green's functions, in the non-autonomous case. The approach we 
shall develop here, by-passes the relative sophistication of the classical approach based 
on evolution families and extends, moreover, to differential-algebraic cases and allows to 
include memory effects, in a simple unified setting. 

To keep the presentation self-contained we construct the Hilbert space setting in sufficient 
detail and formulate our results so that the autonomous case re-appears as a special case 
of the general non-autonomous situation. 

In order to formulate the problem class rigorously and to avoid at the same time to incur 
unnecessary regularity constraints on data and domain for prospective applications, it is 
helpful to introduce suitable extrapolation Hilbert spaces (Sobolev lattices). This will 
be done in the next section. In Section [2] we shall describe a class of non-autonomous 
evolutionary equations and its solution theory. The paper concludes with the discussion 
of a particular application to a class of Kelvin- Voigt-type models for visco-elastic solids to 
exemplify the theoretical findings. 

In the following let H he a Hilbert space with inner product (-I-) and induced norm | ■ |. 

1 Preliminaries 

In this subsection we recall the construction of a short Sobolev chain associated with a 
normal, boundedly invertible operator on some Hilbert space H as it was presented 
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1 Preliminaries 



for instance in |1H Section 2.1]. We note that the construction can be generahzed to the 
case of a closed, boundedly invertible operator on some Hilbert space. However, since the 
operators we are interested in are all normal, we may reduce ourselves to the slightly easier 
case of normal operators. 

Definition. Let if be a Hilbert space and A^: D[N) C ii — )• if be a normal operator with 
G g{N). Then the sesqui-linear form 

{■\-)H,iN):D{N)xD{N)^C 

if,g)^{Nf\Ng), 

defines an inner product on D{N) and due to the closedness of A^, D{N) equipped with 
this inner product becomes a Hilbert space. We denote it by Hi{N). Moreover, by 

{■\-)h.,(n)--HxH^C 

we define an inner product on H and we denote the completion of H with respect to the 
induced norm by H-i{N). The triple {Hi{N), Hq{N), H^i{N)) is called short Sobolev chain 
associated with N, where Ho{N) := H . 

Remark 1.1. Note that the above construction can be done analogously for the /c-th power 
of N for /c G N. The resulting sequence of Hilbert spaces (iifc(A^))^gg , where we set 
Hk{N) := Hi {N'') and H_k{N) := ii_i (A^'^) for A; G N \ {0} is called the Sobolev cham 
associated with N . 

A simple estimate shows that Hi{N) is densely and continuously embedded into Hq{N), 
which itself is densely and continuously embedded into H_i[N). Thus we arrive at a chain 
of dense and continuous embeddings of the form 

iii(Ar)-^iio(Ar)-^ii_i(A^). 

We can establish the operator A^ as a unitary operator acting on this chain. 

Proposition 1.2 (|lli Theorem 2.1.6]). The operator N : Hi(N) -> Hq{N) is unitary. 
Moreover, the operator N : D{N) C Ho{N) H_i{N) has a unitary extension to iio(A^)a 

Since A^ is boundedly invertible its adjoint A^* has a bounded inverse, too. Hence we can do 
the same construction of a Sobolev chain associated with A^*. However, since A^ is assumed 
to be normal we have that D{N) = D{N*), as well as |A^a;| = \N*x\ for each x G D{N). 
This yields that the Sobolev chains for A^ and A^* coincide. Hence, by Proposition 11.21 we 
obtain the following result. 

^We will not distinguish notationally between the operator N and its unitary realizations on the Sobolev 
chain associated with N. 
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Corollary 1.3. The operator N* : D{N*) C Hi:{N) — )• Hk^i{N) has a unitary extension 
to HkiN) for k e {0, 1}. 

Again this realization will be denoted by A^* although this might cause confusion, since 
TV* could be interpreted as the adjoint of the unitary realization of the operator as in 
Proposition 11.21 The adjoint would then be the respective inverse. 

Note that for normal N E L [H, H) =: L{H), the space of continuous linear mappings from 
H to H, we have that 

Hk (N) = H 

as topological linear spaces with merely different inner products (inducing equivalent 
norms) in the different Hilbert spaces Hk{N), k G { — 1,0,1}. This indicates that con- 
sidering continuous linear operators A^ does not lead to interesting chains. 

In the remaining part of this section we consider a particular example of a normal operator 
and its associated Sobolev chain, namely the time derivative in an exponentially weighted 
L^-space (see [HI [121 E] for more details). 



For ^) G M we consider the Hilbert space 

i/,,o(M) := { / e LL(M) I (x ^ exp{-gx)f{x)) G L\R)} 
equipped with the inner product 

{f\9)H,,om '■= j fixTg{x)exp{-2gx) dx {f,g E i/,,o(K)) . 

We define the operator do^g on ifg o(I^) the closure of 

^o,,lc^(M): C'oo(K) C if,,o(M) ^ i/,,o(K) 

where by C'oo(M) we denote the space of arbitrarily often differentiable functions on M with 
compact supporlH. In this way we obtain a normal operator with D\t do^g = q. Hence, for 
7^ the operator 9o,p is boundedly invertible and one can show that ||c}(^^||^|,^ ^^^^ = 

l/|f)|. Thus for ^ 7^ we can construct the Sobolev chain associated with Oq^q and we 
introduce the notation ifg^fc(M) := Hk{do^g) for 7^ and k G {—1,0,1}. For Jmdo^g 
we have as a spectral representation the Fourier-Laplace transform C^: Hg^oiM.) — 
given by the unitary extension of 

Coo (K) C iJ,,o(M) ^ L 



^ (^^^= / exp(-ix|/) exp{-gy)(f){y) dy) 
V v27rjK / 



^The domain of 9o.e consists precisely of the functions / e Hgfi{K.) with distributional derivative lying 
in i/e,o(R)- 
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2 Space-time evolutionary equations 



In other words, we have the unitary equivalence 

where m denotes the multiphcation-by-argument operator in (M) with maximal domain, 
i.e. (m/) (x) := xf{x) for a.e. x G M, / G D{m) := {g G L'^(R) \ {x xg{x)) G ^^(IR)}. 
The latter yields 

do^g = CI (im + g) C^. 
For g ^ we can represent the resolvent do^^ as an integral operator given by 

/X 
u{t) dt {x G Ma.e.) 
-oo 

if g > and 

POO 

(dQ^lu) (x) = - u (t) dt (x G M a.e.) 

J X 

if g < for all u G Hgfi{W). Since we are interested in the (forward) causal situation (see 
Definition 12. 141 below) . we assume g > throughout. Moreover, in the following we shall 
mostly write Oq for Oq^^ if the choice of g is clear from the context. 

Let now denote a normal operator in a Hilbert space H with in its resolvent set. Then 
has a canonical extension to the time- dependent case, i.e., to H^o (M; H) = Hq^o[M.)^H, 
the space of if-valued functions on M, which are square-integrable with respect to the 
exponentially weighted Lebesgue-measure. Analogously we can extend do to an operator 
on o(IR; H) in the canonical way. Then Oq and A^ become commuting normal operators 
and by combining the two chains we obtain a Sobolev lattice in the sense of |1H Sections 
2.2 and 2.3] based on [do, N) yielding a family of Hilbert spaces 

with norms | ■ | ^ ^ ^ given by 

for k,s E { — 1,0, 1}. The operators do and A^ can then be established as unitary map- 
pings from Hi,^k(^;Hs{N)) to H^^k-i C^; HsiN)) for k G {0,1}, s G {-1,0,1} and from 
H,^k (M; H,{Nj) to H,,k (M; Hs-i{N)) for A; G {-1, 0, 1}, s G {0, 1}, respectively. 



2 Space-time evolutionary equations 

Well-posed ness for a class of evolutionary problems 

We are now ready to rigorously approach the well-posedness class we wish to present. We 
shall consider equations of the form 

{doM{mo,d^') + A)u = F, 
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where for simplicity we assume that A is skew-self adjoint in a Hilbert space H and M{t, ■) 
is a material law function in the sense of |10] for almost every t G M. More specifically we 
assume that M is of the formn 



M(mo, ^o"^) = Mo(mo) + do^Mi{mo), 

where Mq, Mi G L(H)), the space of strongly measurable uniformly bounded func- 

tions with values in L{H). We understand 9oMo(mo) + Mi (mo) + A as an (unbounded) 
operator in Hg^Q{M.; H) with maximal domain 

V:={ue H^^o(R; H) \ doMo{mo)u + Au e i/f,,o(K; H)}. 

Note that for u G Hg^Q{W; H) we have 

doMo{mo)u + Mi{mo)u + Au e i/^,,_i(M; H^i{A + 1)) 

and for this to be in H^^q{W; H) is the constraint determining the maximal domain V. 

Hypotheses. We say that T G LfCR; L{H)) satisfies the property 

(a) if T{t) is selfadjoint (t G M), 

(b) if T(t) is non- negative (t G M), 

(c) if the mapping T is Lipschitz-continuous, where we denote the smallest Lipschitz- 
constant of T by |T|Lip, and 

(d) if there exists a set iV C M of measure zero such that for each x & H the function 

R\N 3t^ T{t)x 

is different iabl^. 

Lemma 2.1. Assume that Mq satisfies properties and Then for each t G M the 
mapping 

Mo{t): H ^ H 

■(Mo(-)x)'(t), tGM\iV, 

0, teN 



is a hounded linear operator with ||Mo(t)||L(H) < |Mo|Lip- Thus, Mq G Lf{M.]L{H)) gives 
rise to a multiplication operator MQ{mQ) G L{H^^q{W; H)) given by 



(Mo(mo)«)(t) := (Mo(-)w(t))' (t) 



"^Note that ($ (mo) (t) = ^{t)ip{t) for t e M, e Coo(M;i7) and $: M ^ L{H) strongly measurable 
and bounded. Hence, $ (mo) G L {Hg,o{R; H)) and ||*(mo)||^(^^ ^^(r.^)) < sup^gg ||<I>(t)||^(^) for each 
g>0. 

^If H is separable, then the strong differentiability of T on R \ for some set N of measure zero already 
follows from the Lipschitz-continuity of T. 
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2 Space-time evolutionary equations 



for u G ifp^o(IR; H) and almost every t G M. Moreover, for u G D{do) the product rule 

doMo{mo)u = Mo{mo)u + Mo{mo)dou (2.1) 

holds. In particular, Mo(mo) G L (if^ _i(R; H)) . If, in addition, Mq satisfies property 
then Mo (mo) is self adjoint. 

Proof. Let t ^W\N. The linearity of Mo(t) is obvious. For x G if we estimate 

y^|(Mo(t + /i) -Mo(t))x| < |Mo|Lip|a:| 

for each h e R \ {0}. Thus, \Mo{t)x\ = | (Mo(-)x)' (t) | < |Mo|Lip|a;|, which shows that 
Mo(t) G iv(ii) with ||Mo(t)||L(_H-) < |Mo|Lip- Assuming property (a), we see that the selfad- 
jointness of Mo(t) follows from 

{Mo{t)x\y) = hm ^((Mo(t + h) - Mo{t)) x\y) 

= \iml{x\{Mo{t + h)-Mo{t))y) 
= {x\Moit)y) 

for each x,y E H. It is left to show the product rule (12. ip . To this end, let (p G C*oo(M; H). 
Then we compute 

Mo{t + h)^{t + h)- Mo{t)(f){t) = Moit + h) {4>{t + h) - 4>{t)) + (Mo(t + h) - Mo(t)) 0(t) 
for each t, /i G M. This yields 

]-{Mo{t + h)(P{t + h)-Mo{t)<l){t)) 
h 

= Mo{t + h) (0(t + /i) - 0(t))^ + (^i (Mo(t + /i) - Mo(t))^ (Pit) 

for every t G M, /?- G M \ {0}. The term on the right-hand side in the latter formula tends 
to 

Mo(t)0'(t) + (Mo(-)0(t))' (t) = (Mo(mo)9o0) (t) + (Mo(mo)0)(t) 

for tGM\iV as /i—)-0. Thus, the left-hand side is differentiable almost everywhere and 
since Mo{mo)(p + Mo{mo)do(p G Hgfi(K; H) we derive that Mo(mo)0 G D{do) and 

doMo{mo)(p = Mo(mo)</' + Mo(mo)9o</'. 

The product rule (12. ip for functions in D{do) now follows by approximation. To show that 
the operator Mo(mo) can be established as a bounded operator on Hg^^i(R; H) we observe 
that 
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and thus, 

\M{mo)u\Hg^_,(m;H) = \dQ^ M{mo)u\H,^o{M.;H) 

< \dQ^Moimo)dQ'^u\H^,o(R;H) + \Mo{mo)dQ^u\H^,o(K;H) 

< Q|Mo|Lip + ||Mo(mo)||L(H,,o(IR;H))^ |m|h,,_i(R;/^) 

for each u E H^fi{R; H). □ 

Remark 2.2. Note that the product rule 

aoMo(mo)0 = Mo(mo)0 + Mo(mo)9o0 

for G D{do) can be extended by continuity to G Hgfi{M.] H). Indeed, both the oper- 
ators 9oMo(mo) and Mo(^o) + MQ{mo)dQ are densely defined continuous mappings from 
H^fi{R; H) to i/j, _i(M; H) coinciding on the dense subset i?£,,i(M; H) C iJp,o(K; ^^)- 

Corollary 2.3. Assume that Mq satisfies properties ^ and T/ien 

P = {n G i/,,o(K; i^) I 9oMo(mo)n + An G i/,,o(M; i/) } 
= {n G i/,,o(K; H) I Mo(mo)a*n - An G i/,,o(K; H) } . 

Moreover, we have (Mo(mo)9Q — A)n = (— 9oMo(mo) + 2^Mo(mo) + Mo(''^o) — for all 

uev. 

Proof. Let n G i/g,o(K;//). Then n G P if and only if (aoMo(mo) + A)u G ifg,o(K;i?). 
Since (-2^Mo(mo) - Mo(mo))^^ e ^^e,o(IR; ^^), we have that {doMo^mo) + A)u G Hg^iR] H) 
if and only if 

i7,,o(K; 3 (5oMo(mo) + A)u + (-2f?Mo(mo) - Mo(mo))n 
= Mo(mo)9oM + Au — 2f)Mo(mo)n 
= -(Mo(mo)9* - A)n, 

where we have used that = —do + 2g, which can be verified immediately. □ 

Corollary 2.4. Let e, g > 0. Assume that Mq satisfies the properties and Then 
for u G Hgfl{R] H) we have that 

(1 + edor^doMo{mo)u = doMo{mo){l + edoY^u - edo{l + £ao)-^Mo(mo)(l + edoY\. 

Proof. For u G ifg^o(IK; -f^); "we compute, invoking the product rule (12. ip . that 

((1 + £9o)"^9oMo(mo) - 9oMo(mo)(l + ^^o)"^) n 

= (1 + edoY^ (aoMo(mo)(l + e^o) - (1 + edo)doMo{mo)) (1 + £9o)-^n 

= (1 + edo)-^ idoMo{mo)edo - edo {doMo{mo))) (1 + sdo)-'u 

= (1 + edo)-'edo (Mo(mo)9o - 9oMo(mo)) (1 + edo)-\ 

= edoil + sdor\-Moimo))il + edor'u. □ 
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2 Space-time evolutionary equations 



In the spirit of the solution theory in [TTl Chapter 6], we will require the following positive 
definiteness constraint on the operators Mo,Mo and Mi: there exists a set A^^i C M of 
measure zero with N (1 Ni such that 

3co > 0, ^0 > Vt G M \ iVi, ^ > ^0 : QMo{t) + ^MqW + Mi(t) > cq. (2.2) 
From this we derive the following estimate. 

Lemma 2.5. Assume that Mq satisfies properties (^-|^. Assume that inequality 112. 2\} 
holds and let g > Qq. Then for u G D{do) fl D{A) and a G M we have that 



mt{doMo{mo)u + Mi{mo)u + Au\u){t)e~^^^ dt > Cq / \u{t)\h"^^^ dt. (2.3) 



For the proof of the lemma, we need the following. 

Lemma 2.6. Let g > 0. Assume that Mq satisfies the properties and Then 

for u G D{do) n D{A) and a G M the following equality holds: 



Dlc{doMo{mo)u + Mi{mo)u + Au\u){t)e~^<^* dt 



^(n(a)|Mo(a)n(a))e-2^'"+ / g{u{t)\Mo{t)u{t))e-^^' dt (2.4) 



a 

^Mo{t)u{t) +mcMi{t)u{t) 



{t))e~^''dt. 



Proof. Let u G D{do) fl D{A) and a G M. Note that, since A is skew-selfadjoint, 

^Rc{Au\u){t) = 
for almost every t G M. Hence, the left-hand side in (12. 4p equals 



d\c{doMo{mo)u + Mi{mo)u\u){t)e^'^^^ dt. 
Using the product rule (12.11) and the selfadjointness of Mo(t) for almost every t G M we get 
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that 

2 / 9^e(«9oMo(mo)u + Mi(mo)n|n)(t)e"^^*dt 



-oo 

a 



2 / lHe(Mo(t)9ou(t)|n(t))e~2f* dt + 2 / (Mo(t)n(t) + 91e Mi(t)u(t)|n(t))e-2^* dt 



9=le(aoM(t)|Mo(t)M(t))e"2^* dt 

a 

+ / ^z{Mo{t)dou{t)\u{t))e-'^^' &t + 2 I {Mo{t)u{t) + Uit Mi{t)u{t)\u{t))e-'^'^' dt 



— oo — oo 

a a 



M(-)|(Mo(mo)M)(-))'(t)e-2^*dt- / !He(M(t)| (9oMo(mo)M) (t))e-2^* dt 



+ J mt{Mo{t)dou{t)\u{t))e-'^'' dt + 2 J {Mo{t)u{t) + DKt Mi{t)u{t)\u{t))e-'^^' dt. 

~oo — c 

Using again the product rule (12. ip we obtain 

a 

- I fHe(n(t)| (9oMo(mo)n) (t))e-2^* dt + / ^z{MQ{t)dQu{t)\u{t))e-'^^' dt 



{MQ{t)u{t)\u{t))e^^'' dt. 
Hence, we arrive at 

a 

2 / fHc(aoMo(mo)u + Mi(mo)u|u)(t)e-^^*dt 



— oo 

a a 

■ 1 



n(-)|(Mo(mo)n)(.))'(t)e-2^*dt + 2 / { -Kh{t)u{t) + M^{t)u{t) 



M(a)|Mo(a)M(a))e-2^'^+ / 2Q{u{t)\MQ{t)u{t))e~'^'' dt 



n(t)^e-2f* dt 



a 

2 j (^MQ{t)u{t) + Mi{t)u{t) 



(t) ) e-'^* dt, 
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2 Space-time evolutionary equations 
where we have used integration by parts. 



□ 



Proof of Lemma \2. 51 Using Lemma [2^ and the fact that Mo(a) is non-negative we end up 
with 

a 

^Rt{doMQ{mo)u + Mi{mo)u + Au\u){t)e-^^^ dt 



-oo 

a 



> j ^(^^Mo(t) + ^Mo(t) + fHeMi(t)^ u{t) 



u{t) ) dt 



>co / dt □ 



Our next goal is to show that (12 .Sp also holds for elements in T). For doing so, we need to 
approximate elements in T) by elements in D{do) fl D{A) in a suitable way. 

Lemma 2.7. For each u E Hgfi(R] H) we have that (1 + edQ)^^u — )■ n as e — )■ + . 

Proof. Since the operator family ((1 + £:9o)~^)^>o uniformly bounded, it suffices to note 
that 

(1 + edaY^u - M = (1 + edoY^edQU 
as e 0+ for every u G H^ i (M; H) . □ 

Remark 2.8. It should be noted that literally the same result holds true for do replaced by 
9q. The proof follows with obvious modifications. 

Lemma 2.9. Let e > and let u e V. Then (1 + edoy^u G D{do) fl D{A) and the 
following formula holds 

(1 + edo)~' {doMo (mo) + Mi (mo) + A)u 

= {doMo (mo) + Ml (mo) + A) {1 + e^o)"' u - edo (1 + edo)'^ Mo (^o) (1 + edo)'^ u 
+ (1 + edoY^ Ml (mo) u - Mi{mo) (1 + edo)'^ u. (2.5) 

Moreover, we have 

{doMo (mo) + Ml (mo) + A) {1 + edo)'^ u -r [SoMq (mo) + Mi (mo) + A)u {e ^ 0+) 
in F^,o(K; H). 
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Proof. With the help of Corollary 12.41 and the fact that A and (1 + eOq) ^ commute, the 
formula f l2.5p follows. From (12. 5p we read off that (1 + edo)~^u G V. Moreover, since 
(1 + edoy^ u e Hg^i (M; H) we have doMo (mo) (1 + edo)'^ u e H^^o (M; H). Defining 

(doMo (mo) + Ml (mo) + A) {1 + e^o)"' u =: F e H,^o (M; H) , 

we get that 

F - {doMo (mo) + Ml (mo)) (1 + edo)-' u = A{l + edo)-' u e i/,,o (M; H) . 

Hence, {1 + edo)~^ u G D (A) . According to Lemma 12.71 the left-hand side of equation 
(12. 5p converges to {OoMq (mo) + Mi (mo) + A) u and the last two terms on the right hand 

side cancel out as e — )■ + . Moreover, since (edo (1 + edo)~^ Mq (^o) (1 + £do)~^ u) 

\ / £>0 

is bounded in Hgfi{M.] H), there exists a weakly convergent subsequence. Using that 
edo (1 + edo)~^ Mq (mo) (1 + edo)~^ m — )■ in Hg_i{R] H), we deduce that 

edo (1 + edoy^ Mo (^o) (1 + edoy^ u^O {e ^ 0+), 

and thus 

(^oMo (mo) + Ml (mo) + A) {1 + edo)~^ u {doMo (mo) + Mi (mo) + A)u (e -> 0+). 

□ 



The most important step to generalize the statement of Lemma [231 to the case of elements 
in T> is the following result. 

Lemma 2.10. Assume that Mq satisfies the properties and Let g > 0, a E MUjoo} 
and let G : Hg^Q{W; H) — )■ M 6e continuous. Moreover, assume that for u G D{do) fl D{A) 
we have that 

a 

J D\c{{doMo{mo) + Mi{mo) + A) u\u){t)e-^<'' dt > G{u). 

—oo 

Then the latter inequality holds for all u eT). 

Proof. Let u eT). According to Lemma [2.91 we have that 

(^]-oo,a](?Tlo)M|(5oMo(mo) + Ml (mo) + A)u) Hg,o{^-H) 

= lim(X]_oo,a](?^o)(l + £:9o)"^M|((9oMo(mo) + Mi (mo) + A){1 + £:9o)"^M)Hg,o(iR;H), 

where we have used that the multiplication operator X]^oo,a]{fno) with the cut-off function 
^]-oo,a] is a bounded operator on Hgfi{M.;H) and that (1 + edo)^^ converges strongly to 
1, by Lemma 12.71 With the assumed inequality and the fact that for e > we have 
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2 Space-time evolutionary equations 

+ e D{do)nD{A) = H^^i(R; H) n H^^o(R; Hi{A + l)), by Lemma[231 we obtain 

that 

^e{X]-^oo,a]{'mo)u\{doMo{mo) + Mi (mo) + A)u)h^^o{K;H) 

= lim9ic(X]_oo,a](rno)(l + edoy'u\{doMo{mo) + Mi (mo) + A){1 + edo)-\)H,,,(R;H) 
>limG((l + edoy\) = G(u). □ 

e-->-0 

Corollary 2.11. Assume that Mq satisfies properties (^-([^. Assume that inequality ^2.2\] 
holds and let g > Qq. Then for u eT) and aGMUjoo} we have that 

a a 

0^e((9oMo(mo) + Mi(mo) + A)M|M)(t)e-2^*dt > Co ! \u{t)\^e-'^'' &t. 



Proof. The statement is immediate from the Lemmas 12.51 and 12.101 with 



a 



G{u)=Cq I |u(t)pe ^^*dt = Co(X]_oo,a]("^o)^i|^i)Hg,o(M;/f)- □ 



-oo 



Lemma 2.12. Assume that Mq satisfies the properties (^-([^j and that inequality (dj 
holds. Let g > go and u E {v e Hg^iR; H) \ {Mo{mQ)d*Q - A)v E H^fi{R; H)}. Then the 
inequality 

j 'nt{{Mo{mo)d; + Mi{moy-A)u\u){t)e-^''dt>Coj \u{t)\'^e-^'' dt 

R R 

holds. 



Proof. By Corollary I2.3[ we deduce that u E T) and that for a G M we have 
mt{{Mo{mo)d*o + MiimoY - A)u\u){t)e~'^'' dt 

a 

I fHc((-aoMo(mo) + 2^Mo(mo) + Mo(mo) + Mi (mo)* - A)u\u){t)e-^s' dt. (2.6) 
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With Lemma we get for u e D{do) n D{A) and a G M that 

mt{{-doMo{mo) + 2^Mo(mo) + Mo(mo) + Mi (mo)* - A)u\u){t)e-^'' dt 



--{u{a)\Mo{a)u{a))e~^^'' - / g{u{t)\Mo{t)u{t))e~^'' dt 



u{t) ) e"'^* dt 



+ / (2^Mo(t)n(t) + Mo{t)u{t)\u{t))e^'^'' dt 



— {u{a)\Mo{a)u{a))e-'^<^'' + / Q{u{t)\MQ{t)u{t))e-^^' dt 



a 

+ j (^MQ{t)u{t)+^tMi{t)u{t) u{t)'^e~^^' dt 



> --(n(a)|Mo(a)n(a))e-2^"^ + co / {u{t)\uit))e-^^' dt 



Letting a — j- oo, we deduce that 



$Ke((-9oMo(mo) + 2gMo{mo) + Mo(mo) + Mi (mo)* - A)u\u){t)e-^^' dt 



> Co / (u(t)|M(t))e-2^* dt. 



Now, Lemma [2.101 imphes the latter inequahty to hold for all u eV. The assertion follows 
from equation (12. 6p . □ 

Theorem 2.13 (Solution Theory). Let A: D{A) ^ H ^ H be skew- self adjoint and 
Mo, Ml G L^(]R; L(i7)). Furthermore, assume that Mq satisfies the hypotheses (^-|^ and 
that 1(2. 2\] holds. Then the operator OqMq (mo) + Mi (mo) + A is continuously invertible in 
Hg ()(M', H) for each g> Qq. A norm bound for the inverse is 1/cq. Moreover, we get that 

{doMo (mo) + Ml (mo) + A)* = (Mo (mo) 9* + Mi (mo)* - A) , 

where the latter operator is considered in Hg^Q{M.; H) with maximal domain. 
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2 Space-time evolutionary equations 



Proof. Let g > Qq. By Corollary 12 . 1 1 1 we have that 

fRt{u\{doMo{mo) + Mi(mo) + A)u)Hg,o{R;H) > Co{u\u)Hg,oiR;H) 

for all u E V. This implies that the operator d^Mo (mo) + Mi (mo) + A has a continuous 
inverse with operator norm less than or equal to the constant I/cq. 

It remains to show that 9oMo(mo) + Mi (mo) + A maps onto Hgfi{R;H). For this we 
compute the adjoint of 

B := (doMo (mo) + Mi (mo) + A) 

considered as an operator in Hgfi{R; H). Let f e D{B*) C H^^o(R; H). Then for all u G P 
and e > we obtain with the help of equation fl2.5p that 

{Bu\{{l + edorrf)H,,i^;H) 

= ((l+.9o)-^i3n|/>^^^^^(^^^) 

= {B{1+ edo)-\\ />^^ ,^(iR.^) + {-edo (1 + edo)"^ Mq (mo) (1 + edo)"^ u\f)H,M^;H) 
+ ( (1 + edo)-^ Ml (mo) u - Mi (mo) (1 + edo)~^ u\ />^^_^(k.h) 



(2.7) 



( (1 + edor u\ ^7)^^ + ( - (1 + edor Mo (mo) (1 + ^^o)"' )7 



+ (u 



((1 + edo)-' Ml (mo) - Mi (mo) (1 + e^o)"')* / 



Hence, we deduce that (1 + ed^)-^[D{B*)] C D{B*) for e > 0. Moreover, we have (1 + 
ed^y^f e D{d^). Thus, for u e H^^i{R; Hi{A + 1)) C D{B) and £ > we get that 

<«|i3*(l + .9*)-V>^^,„(,^^) 

= { {doMoimo) + Ml (mo) + A) u\ (1 + ^5*)"V>^^ 

= {u\{Mo{mo)dl + Ml (mo)*) (1 + ^9o*)"V>^^,.,(^^^) + (^^(1 + ^^o)'V >^^,(^^^) • 

Since i^'e,i(M; Hi{A + 1)) is a core for A, we deduce that (1 + ed^y^f E D{A) for e > 0. 
Moreover, we have 

B* (1 + 6d*r' f = (Mo(mo)9* + Ml (mo)* - A) (1 + ed*)"' f. 

Using (12. 7p we can estimate 



B*il + ed*r'f 



Ih„o(M;H) =SUp{|('«|^*(l+^5o) V>H„o(M;H) 
+ 2\\Mi{mo)\\L{Hg,o{R;H))\f\Hg,oiR;H) 



u eV, \u\ 



< 1 
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for every e > and thus, we find a weakly convergent subsequence in Hg^Q(R;H) C 
ifg,_i(M; i/)ni/g,o(K; H^i{A+l)). Moreover, note that {Mo{mo)d^ + Mi (mo)* - A) (1 + ed^) 
converges to {Mo{mo)d^ + Mi (mo)* - A) f in H^^^i{R; H)r]H^fi{R; H^i{A + l)). Thus, by 
the (weak) closedness of B* we derive 

B*f = (Mo(mo)9o* + Ml (mo)* - A) f 

and 

D{B*) C{fe iJ,,o(K; H) \ (Mo(mo)9* + Mi (mo)* - A) f e i/,,o(K; H) } . 
We define 

C:D{C)CH^^o{^;H)^H,^o{^;H) 

f ^ (Mo(mo)9* + Ml (mo)* - A) /, 

where D{C) := {/ G i/£,,o(M; if) | (Mo(mo)ao - A) f e Hg^oi^; H)} . Lemma nA2\ ensures 
that C is one-to-one. Thus, so is B*. According to the projection theorem we have the 
orthogonal decomposition 

if,,o(M; H) =N{{doMo{mo) + Mi (mo) + A)*) © i?(9oMo(mo) + Mi (mo) + A) 
={0} © R{doMo{mo) + Ml (mo) + A) 

and this establishes the onto-property of c}oMo(mo) + Mi (mo) -|- A. Moreover, we get that 
(i5-i)* = (i3*)-i c C-\ The first operator is left-total. Thus, B* = C. □ 



Causality 

At first we give the definition of causality in our framework. 

Definition 2.14. Let if be a Hilbert space, g > and G: D{G) C Hgfi{R;H) 
Hgfl(R; H). Then G is called (forward) causal, if for each a G M and each f,g& D{G) the 
implication 

X]„oo,a](mo)(/ - (?) = ^ X]_oo,a](mo) {G{f) - G{g)) = 

holds. 

Now, we want to show that our solution operator (9oMo(mo) + Mi (mo) + A)~^ is causal 
in ii,,o(M; H). 

Theorem 2.15 (causal solution operator). Under the assumptions of Theorem \2.1'A the 
solution operator (9oMo(mo) + Mi (mo) + A)~^ is causal in Hgfl{R] H) for each g > g^. 



19 



2 Space-time evolutionary equations 



Proof. Let / G //j,,o(K; H) with X]_oo,a]("^o)/ = 0. We define u := ((9oMo(mo) + Mi (mo) + 
y4)~^/ G P and estimate, using Corollary 12.1 ![ 



9^e(/|M)(t)e"2f* dt > Co / \u{t)\^e~^^' dt, 



which shows that X]_oo,a]("^o)^ = 0. Thus, due to linearity, the solution operator is causal. 

□ 



An illustrative example 

To exemplify what has been achieved so far, let us consider a somewhat contrived and 
simplistic example. 

The starting point of our presentation is the (1 + l)-dimensional wave equation 

OqU - dfu = f onRxR. 
As usual we rewrite this equation as a first order system of the form 

In this case we can compute the solution by Duhamel's principle in terms of the unitary 
group generated by the skew- self adjoint operator 

-^i 
-di 

This would be the simplest autonomous case. Let us now, based on this, consider a slightly 
more complicated situation, which is, however, still autonomous: 



u 

XM\]_e,£[(mi) y ' \^ X]_£,e[(mi) y ' \-di ) ) \v 



do'f 




(2.9) 



where X/(mi) denotes the spatial multiplication operator with the cut-off function X/, i.e. 
(X/(mi)/) {t,x) = Xi{x)f{t,x) for almost every (t,x) G R x M, every / G H^fi{R;L\R)) 
and / CM. In the notation of the previous section we have 



Mo (mo) := 



^R\]-e,0[("^l) 

XR\]-e,el{mi] 
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and 



and both are obviously not time-dependent. Note that our solution condition (12. 2p is 
satisfied and hence, according to our findings, problem ( 12. 9 p is well-posed in the sense of 
Theorem 12.131 By the dependence of the operators Mo(mo) and Mi (mo) on the spatial 
parameter, we see that (12. 9p changes its type from hyperbolic to elliptic to parabolic and 
back to hyperbolic and so standard semigroup techniques are not at hand to solve the 
equation. Indeed, in the subregion ] — e, 0[ the problem reads as 



u 



-OA f u\ _ fd^'f 



V J \ —di J \v J \ 
which may be rewritten as an elliptic equation for u of the form 

u - d\u = Oq^J. 

For the region ]0, £:[ we get 



which yields a parabolic equation for u of the form 

dou - dlu = V- 

In the remaining subdomain M\] — e,e[ the problem is of the original form (12. Sp . which 
corresponds to a hyperbolic problem for u. 

To turn this into a genuinely time-dependent problem we now make a modification to 
problem (12. 9p . We define the function 

{0 if t < 0, 
t ifO<t<l, (teR) 
1 if 1< t 

and consider the material-law operator 



Mo(mo) = v^(mo)f'^^\l-°[("^^ 





XK\]_e,£[(mi) 

which now also degenerates in time. Moreover we modify Mi (mo) by adding a time- 
dependence of the form 



Ml (mo) 



^]-oo,0[("^o) + ^[O,oo[("^o)^]-e,0[("^l) 

X]_oo,o[("^o) + ^[o,oo[("^o)^]-£,£[("^i) 
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2 Space-time evolutionary equations 



We show that this time- dependent material law still satisfies our solvability condition. To 
this end let g > 0. Note that 



and thus, for t < we have 



iftG]0,l[, 
otherwise 



pMo(t) + ^Mo(t)+5HeMi(t)= (j J ) > 1- 

For < t < 1 we estimate 

gMo{t) + ^Moit)+mcM^{t) 



2 " ; V XR\]_,,,[(mi) J ' \ X]_,,,[(mi) J - 2 

and, finally, for t > 1 we obtain that 
gMo{t) + ^Mo{t)+<yicMi{t) 

Remark 2.16. Note that the spatial operator ( ^ J in the previous example can be 



-di 

substituted by every skew-selfadjoint operator. In applications, it turns out that this oper- 

/ Q C*\ 

ator typically is a block operator matrix of the form I ^ n I ' where C is a densely de- 



^ -C ^ 

fined closed linear operator between two Hilbert spaces. Indeed, even the one-dimensional 
transport equation shares this form, if one decomposes the functions in their even and odd 
parts (see [T5| p. 17 f.]). Moreover, it should be noted that the block structures of the 
/ g C*\ 

operator ( (j q j ^^'^ operators Mo(mo) and Mi (mo) need not to be comparable 

(it turns out that this naturally arises in the study of boundary control systems, cf. |14] . 
|15|). In those cases the semi-group approach for showing well-posedness is not applicable, 
without further requirements on the block structures of the involved operators. 



Some perturbation results 

In applications, it is useful to have a perturbation result at hand. To this end, we assume 
we are given a linear mapping 

Moo : D{M^) C Pi H,,o (M; H) ^ f] iJ,,o (M; H) 
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for some > in the way that for all > f?o we have that D{M^) C HgQ{M.;H) is 
dens^ and Moo considered as a mapping from Hgfi(R; H) to Hgxi(R; H) is continuous. 
The assumptions give rise to a continuous extension, denoted with the same symbol. A 
straightforward consequence of our previous findings is the following. 

Theorem 2.17. Let A: D{A) C H H be skew- self adjoint and Mo,Mi E Lf{R; L{H)). 
Furthermore, assume that Mq satisfies the properties and that holds. Assume 

that 

limSUp||Moo||i(^^,,(K;iy)) < Co. 

Then there exists Qi > such that the operator doMo (ttiq) + Mi (mo) + M^o + A is contin- 
uously invertible in HgQ(R;H) for each g > Qi- If, in addition, is causal, then so is 
{doMo (mo) + Ml (mo) + + Ay\ 

Proof. Let Qi > he such that ||Moo||j;^(j:^ a{R-H)) ^ '^o for all g > Qi. Let / G Hg^iR; H). 
Then, the mapping 

u ^ (9oMo(mo) + Ml (mo) + A)"' (/ - M^u) 

is a strict contraction, by Theorem 12.131 Observing that u G Hg^Q(R; H) satisfies 

(^oMo (mo) + Ml (mo) + M^ + A)u = f 

if and only if it is a fixed point of $, we get existence and uniqueness of a solution with the 
help of the contraction mapping principle. If Moo is causal, then so is $ as a composition 
and a sum of causal mappings. Hence, (c?oMo (mo) + Mi (mo) + Moo + ^)~^ is causal. □ 

Remark 2.18. Note that this perturbation result applies similarly to the case of non-linear 
perturbations if the best Lipschitz constant jMool^^p of the perturbation Moo considered 
as an operator in (M; if) , ^) G ]0, cx)[, satisfies 

limsup I Moo I ^ Lip < Co. 

It is possible to derive the following more sophisticated perturbation result, which needs 
little more effort. We introduce the following notation: For a closed subspace V H we 
denote hy cy '■ V H the canonical embedding of V into H . It turns out that then the 
adjoint Ly : H V is the orthogonal projection onto V. Consequently Py '■= i-v^v ■ 
becomes the orthogonal projector on V and 1 — Py = Py± = iyxi^^. 

^Note that as an example CooO^',H) is dense in Hg,o{M.; H) for all g> {). 
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2 Space-time evolutionary equations 



Theorem 2.19. Let A: D{A) C H ^ H be skew- self adjoint and Mo,Mi G L^(R;L{H)). 
Furthermore, assume that Mq satisfies the properties (^-|^. Moreover, assume t t— )■ 
N[Mq (t)) to be time-independent, i.e., for a// 1 G M we have 

N{Mo {t)) = N{Mo (0)) =: V. 

We further assume that for some set of measure zero Ni the following estimates hold: 

3co > 0, ^0 > OVt G M\ A^i : L*yD\t Mi{t)Lv > Cq, (2.10) 

and 

3ci > Vt G M \ iVi : L*y^Moit)Lv± > c^. (2.11) 

Furthermore, assume that limsupg_^oo ll^^oo||i(j:/ o{r-h)) ^^'^ there exist g,e > such 

that for all g > g and u G D{Moo) 

(9^cMooPyM| i^yM)H,,o(R;H) > " ^o) |Pfm|h,,o(R;H) ' 

Then there exists > such that the operator OqMq (mo) + Mi (mo) + Moo + A is contin- 
uously invertible in Hgfi{W;H) for every g> gi. If, in addition. Moo is causal, then so is 
(doMo (mo) + Ml (mo) + Moo + A)-\ 

The resuh foUows by adapting the method of proof of Theorem 12. 131 The crucial estimate 
to conclude the proof of Theorem 12.191 is given in the following lemma. 

Lemma 2.20. Let g > go- Assume that Mq satisfies the properties 1^-^, that t H- 
7V(Mo(t)) is time-independent and that the inequalities ^2.10\) - [2.11\) hold. Then for all 
e g]0, Co[ there exists c > such that for all sufficiently large g and for u ^ V and a G M 
we have that 



((aoMo(mo) + Mi (mo) + A) u\ u) {t)e-^'' dt 

a a 

>QC j |Pyxn(t)|^e-2^*dt+ (co -e) j |Pyu(t)|^ e"^^* dt. (2.12) 

— oo — oo 

Proof. In order to prove fl2.12p observe that by Lemma 12.101 it suffices to verify the in- 
equality for u G D{do) n D{A). Moreover, by Lemma [2. 6 [ we only need to estimate 



a 

^{u{a)\Mo{a)u{a))e-^''' + j (^gMo{t)u{t) + ^Mo{t)u{t) + ^Kt Mi{t)u{t) 

— oo 

Since Mo(a) is non-negative, we are reduced to showing an estimate for 

^Mo(t)0 + ^Mo(t)0 + Mi(t)0 



uit) } dt. 
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for alH e M and (pe H. Using that Mo(t) vanishes onV = N{Mo{0)), we get that 



gMo{t)(f)+^Mo{t)(f) + D\zMi{t)4> 



+ 2 ($HeMi(t)PyX0| Py0) + {mtMi{t)Py±(P\ Py±(P) 
+ {D\zM^{t)Py4>\Py(P) 

> (^QCi - ^|Mo|Lip||Mi(mo)||L(/f„.o(M;/f))^ \Pv^(pf 
- (iMolLip + 2||Mi(mo)|U(H,.o(R;//))) \Pv^<P\ \Pv<l)\ + Co \Pv 



The assertion follows now by applying the trivial inequality 2ab < -^a + 6b for a,b,6 > 



0. 



□ 



Proof of Theorem \2.19{ . We denote B := 9oMo(mo) + Mi (mo) + A + Moo considered as an 
operator in Hgfi{M.;H). Note that the maximal domain in Hq^q(R; H) coincides with V. 
Moreover, since Mqo is continuous, we have, by Theorem I2.13[ that 

B* = Mo{mo)d* + Ml (mo)* -A + M*^ 

with domain being equal to V, by Lemma 12. 3[ Let e > and choose g such that 
supg>^||Moo||L(//,,o{R;H)) < oo and 

{D\tMooPvu\ Pv'«)h,,o(R;H) > - Co) l^y^lH,,o(K;H) 

for each g > ^and u G D{Moo)- Note that due to continuity the latter inequality also holds 
for every u G Hgfi{M.;H). For m G P there exists, according to Lemma l2.20[ a constant 
c > such that for all sufficiently large g >g we have 

DKt{Bu\u)H,,o{^-H) = ((9oMo(mo) + Mi (mo) + A)u\ u)^^^^^^^.^^ + ^Rt {Moou\ u) ^^^^^^^^.^^ 

> ^c|Pyxn|^^^^(K;H) + (cO - I) l^y«lH,,o(M;H) + " Co) l^y«l^^,o{M;H) 

- 2||Moo||L(Hp,o(R;//))|-Pv'"|He,o(IR;J^)|-Pv-LM|//g_o{K;H) 

- ll^oo||L(//p,o(K;//))|-PyiW|/fg_o(M;H) 
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3 An application to a Kelvin-Voigt-type model in visco-elasticity 



for each < 5 < |. By possibly increasing g such that 

1 2 

gC- \\Moo\\L{Hg,oOR;H)) " -^\\Moo\\l(H ,,oi^;H)) > C > 

we deduce that for all n G the estimate 

^Kt{Bu\u)Hg^o{R;H) > c|u|^^_^(K.^) (2.13) 

holds for all sufficiently large g. By D\t{Bu\u)Hg,oi'^-H) = ^ic {B*u\u) jj g(^.H) ^i^) = 
T? = D{B*), we deduce that B* is one-to-one. Hence, B is continuously invertible and onto. 
For showing the causality of B~^ in case of a causal operator Moo it suffices to prove, that 
an inequality of the form 

a a 

«./(H4.>W.-d«>j/M*)|V-d* 

— oo — oo 

holds for every a G M and u eT). The latter can be shown as above, observing that due to 
the causality of M^o we have 

9^c(MooM|X]_oo,a] ("2o)m)h,,o(IR;H) = 5^e(X]_oo,a] {^0) M ^u\Xy ^ ^a] ("2o)m) H,,o(IR;H) 

= 5^e(MooX]_oo,a](?Tio)^^|X]_oo,a](?Tio)n)Hgo(IR;H)- □ 



3 An application to a Kelvin-Voigt-type model in 
visco-elasticity 

Although, applications are obviously abundant by simply extending well-known autonomous 
problems to the time- dependent coefficient case, we intend to give a more explicit applica- 
tion here to illustrate some of the issues that may appear in the non-autonomous case. A 
more straight-forward application would be for example solving Maxwell's equations in the 
presence of a moving body, which reduces via suitable transformations to solving Maxwell's 
equations with the body at rest but coefficients depending on time, [3j . Applying the above 
theory to this case avoids the intricacies of Kato's method of evolution systems employed 
in [3]. As a by-product, the assumptions needed are considerably less restrictive. 

As a more intricate application we would like to elaborate on here, we consider a time- 
dependent Kelvin- Voigt material in visco-elasticity. In [2] such a material is considered 
in connection with modeling a solidifying visco-elastic composite material and discussing 
homogenization issues. We shall use this as a motivation to analyze well-posedness in the 
presence of such a material under less restrictive assumptions. 

In this model we have the equation 

dor]{mo)doU - DivT = / 



26 



linking stress tensor field T with the displacement vector field m, accompanied by a material 
relation of the form 

T={C{m^) + D{m^)do)S (3.1) 

where Div is the restriction of the tensorial divergence operator div to symmetric tensors 
of order 2 and £ := Grad u with Grad denoting the symmetric part of the Jacobian matrix 
d ® u oi the displacement vector field u. The operators C, D and rj are thought of as 
material dependent parameters. Here the case D (mo) = would correspond to purely 
elastic behavior. Introducing v := df)U as a new unknown we arrive, by differentiating 
(IM]), at 

dof]{mQ)v — DivT = / 
do {C (mo) + D (mo) ^o)"^ T = Gradu, 

where we can choose g large enough, such that 

C (mo) + D (mo) do = {C{mo)d^' + D{mo)) do 



gets boundedly invertible. Assuming for sake of definiteness vanishing of the displacement 
u on the boundary as a boundary condition we obtain an evolutionary equation of the form 

^°(''V\c(mo) + D(mo)9o)-0^(-Grad ))(r) = (o)' ^^'^^ 



where the choice of boundary condition amounts to replacing Grad by the closure Grad 
of the restriction of Grad to vector fields with smooth components vanishing outside of a 

2 

3x3,sym 

of {ny © {ny''^ with its natural norm, where the second block-component space 

2 

3x3, sym 



compact subset of Q. The underlying Hilbert space is the subspace {Q)^ © -^1x3 sym (^) 
if {ny © {ny''^ with its natural norm, wl 

-^3x3 sym (^) dcuotcs the restriction of L"^ (Q)"^^^ to symmetric matrices with entries in 
(Q) . Note that then ^ ^ ^ ' ^ is skew-self adjoint (see e.g. [TT| Section 5.5.1]). 

The operator families (C (t))t6M and {D{t))t£w are assumed to be uniformly bounded in 
-^3x3 sym (^) • Further constraint will of course be required to satisfy the assumptions of our 
solution theory above. We are led to a material law operator of the form 

A1(9„-') = ("'"»' ° 



(C(m„) + D(m„)a„)-' 



To deal with the term {C (mo) + D (mo) do)~ we need a projection technique. For this we 
recall that for a closed subspace V of the underlying Hilbert space -^1^3 sym denotes 
the canonical injection of V into L\^^^^^{Vt). The solution theory for the Kelvin- Voigt- 
type model is then summarized in the following theorem. 
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3 An application to a Kelvin-Voigt-type model in visco-elasticity 



Theorem 3.1. Let Q C M.^ be open and V be a closed subspace of L\^^^^^{VL). Let 
C e L^(M;L(L2^3_,y^(fi))), r/ G L^(M; L(L2 (fi)^)) and B G L^(M;L(V^)). Assume that 
0,7] satisfy the properties 1^-^. We set 

for a// 1 G M and we assume the existence of c > such that for all t we have 

mcB{t) > c, L*y^C{t)Lv± > c, ri{t) > c. 

Then for all sufficiently large g we have that for all F G Hq^q{W;L'^ (Q)^ © -^1x3 sym (^)) 
there exists a unique solution {v,T) G H^^qIM.] L"^ {Q)^ © Lf^^^^^ {Q)) of the equation 



'f]{mo) \ f ^ — Div\\ / V 



{C (mo) + D (mo) doy^ J ' V-Grad J J \T 

The solution depends continuously on the data. The solution operator, mapping any right- 
hand side F to the corresponding solution of the latter equation, is causal. 

Proof. The proof rests on the perturbation resuh Theorem 12. 171 Since the top left corner 
in the system under consideration clearly satisfies the solvability condition (12. 2p . we only 
have to discuss the lower right corner. For this we have to find a more explicit expression 
for {C (mo) + D (mo) Oq)'^ . An easy computation shows that 

{C (mo) + D (mo) do)-' 

4C(mo)iy + B{mo)do 4C(mo)iy± 

{L*y^C{mo)Ly±) h*y^C{mo)iv ^ J ^ ^ 

with 

W{mo) 

{LyC{mo)Lv + B{mQ)do ~ iyCimo)Ly±{L'l^^C{'mo)Ly±)~^Ll,^C{rno)iv) ^ 

{Ly^C{mo)tvy^^ 

Denoting ^(mo) := [_^,*^^c{mo)tylr'^v.C{rnoW l) ' 

{C (mo) + D (mo) do)'' = S {mo)W (mo) S {mo)\ 

by the selfadjointness of C(mo). Now, the top left corner of W{mQ) may be expressed with 
the help of a Neumann expansion in the following way 

{B{mo)do + LyC{mo)Lv - 6yC(mo)iyx(ty±C(mo)6yx)~^6yxC(mo)iy) ^ 

= dQ^B{mo)~^ (l + d^^B{mo)~' (4C(mo)i\/ - iyC(mo)tyx(iy±C(mo)iy-L)"^''yxC(mo)ty))~^ 
= d,'B{mor' + d^'M^ 
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for some suitable Mqo, satisfying \\Moo\\ L{Hg o{R;V(bv^)) — >■ as — )■ oo. Thus, we arrive at 
{C{mo) + D{mo)dor' 

+ ^(^o)f^°''^°°!!VM*. (3.3) 



0^ 

The first term on the right-hand side can be computed as follows: 

(o {L^y^Zw^r') ^^"^'^^ = (o (.^xC(mVyx)-i 

For the second and third term on the right-hand side we observe that S G L^(R; L{V © 
V-^)). Moreover, S is Lipschitz-continuous. Indeed, observing that 

< \\{t*y^C{t)Ly.r'\\L^y.^\\t*y4C{s) - C{t))Ly4LiV^)\\{L*y^C{s)Ly.rXiV-) 

<c-'\C\uAs-t\ 
and that, using |C|oo := supjgjg ||C(t)||, 

UL*y^C{t)Ly^)-h*y^C{t)Ly - {i*y^C {s)iy^)-h*y^C {s)Ly\\LiV,V^) 

< ||(.;.xC(t).vx)"iL(yX)||44C(t) -C(s)).v'|U(V^,VX) 

+ \\{L*y^C{t)Ly.r' - (^^xC(s)6yx)-l|U(;.X)||6;.,C(s).y|U(y,V'X) 

< {c-'\C\up + c-^\C\up\CQ\s-t\, 

for s,t G M, we derive the Lipschitz-continuity of S. Thus, S satisfies the hypothesis ([d]), 
since L"^^^ ^^^^{fl) is separable. Hence, we can compute the second term on the right-hand 
side of fl3.3p as follows by using the product rule 02.1 p 

S(m„)(*"^^'"°)"°)s(m„)- 



9o doS{mo) [ " 



5(mo)* 



Defining Mo(mo) := ( ° , ^ ° . Mi (mo) := S{mo) f ^("^o)^' SimoY 



{i*y^C{mo)Ly±y' J' ■ ^ ojy Q Q 

and Moo := S{mo) (^^o'B{mo)-' 0^ ^ ^^^^^ ^ao^^Mo, 0^ ^^^^^.^ ^^^^ 

^0 (C (mo) + D (mo) So)"' = 5oMo(mo) + Mi (mo) + Moo- 
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Clearly, M^o can be considered as a perturbation as in Theorem I2.17[ since 

||^oo||L(//,,o(M;y®y^)) ^0 {g-^ oo). 

Moreover, one easily obtains that y{t{Mi{t)Lv(j)\iv4')v®v^ ^ for all t G M 

and G Since Mo{t) is strictly positive on V-^ and Mo(t) is uniformly bounded in t, 
we get estimate (12. 2p for sufficiently large Qq. Theorem 12.171 concludes the proof. □ 

Remark 3.2. The assumption on the subspace V in the above theorem expresses the fact 
that the null space of D (mo) is non-trivial due to degeneracies in various regions in Q, but 
it is implied that these regions do not vary in time. In other words there may be stationary 
areas in which the material exhibits purely elastic behavior and others showing different 
forms of visco-elastic behavior. 
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